Self-organized Vortex State in Two-dimensional Dictyostelium Dynamics 
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We present results of experiments on the dynamics of Dictyostelium discoideum in a novel set- 
up which constrains cell motion to a plane. After aggregation, the amoebae collect into round 
"pancake" structures in which the cells rotate around the center of the pancake. This vortex state 
persists for many hours and we have explicitly verified that the motion is not due to rotating waves 
of cAMP. To provide an alternative mechanism for the self-organization of the Dictyostelium cells, 
we have developed a new model of the dynamics of self-propelled deformable objects. In this model, 
we show that cohesive energy between the cells, together with a coupling between the self-generated 
propulsive force and the cell's configuration produces a self-organized vortex state. The angular 
velocity profiles of the experiment and of the model are in good quantitative agreement. The 
mechanism for self-organization reported here can possibly explain similar vortex states in other 
biological systems. 



Spontaneous organization of self-propelled particles 
can be found in a variety of systems. Examples in- 
clude the flocking of birds [jlj, the movement of traffic 
\2\ and pedestrians || and the collective motion of ants 
Not surprisingly, the physics of self-propelled par- 
ticles has recently attracted considerable attention ||. 
The systems have been theoretically analyzed using con- 
tinuum, Navier-Stokes like, equations and discrete nu- 
merical models. The numerical models have in common 
that the objects are treated as point-particles. However, 
in a variety of systems the actual shape and plasticity of 
particles play an important role and hence they can not 
be modeled in this oversimplified manner. In this Letter 
we present experimental data in one such system, Dic- 
tyostelium discoideum cells, and introduce a model that 
treats cells as deformable objects rather than point par- 
ticles. Our results provide evidence for a localized vortex 
state in this biological system. 

The developmental dynamics whereby Dictyostelium 
discoideum is transformed from a solitary amoeba state 
to a functional multicellular organism is of interest to 
both biologists and non-equilibrium physicists |J. Past 
efforts have elucidated the nonlinear chemical wave sig- 
nals used to guide aggregation |j]-|9| as well as the chemo- 
tactic instability |l(J which causes a density collapse to 
one-dimensional "streams" of incoming cells. However, 
much less is known regarding subsequent events, espe- 
cially with regard to organized cell motion in the later, 
multicellular stages of the day-long developmental pro- 
cess. 

Dictyostelium cells are grown in liquid media and 
plated without nutrients onto a glass surface [ pT[ . An 
additional thin layer of agarose is then overlaid on the 



cells. This has the effect of restricting cell motion to the 
plane; in fact, the multicellular states that form are at 
most a few monolayers deep. Cells aggregate normally 
and form round "pancake" structures. In Fig. la, we 
have presented a typical snapshot of the system. In each 
of the observed structures, the cells have organized their 
motion into a coordinated vortex state in which they ro- 
tate around the center of the pancake. The rotation can 
be either clockwise or anti-clockwise depending on ini- 
tial conditions and can persist for tens of hours. Fig 
lb shows a close-up of one structure, where now the cells 
have been illuminated by using a strain in which the gene 
for green fluorescent protein jl2j has been fused to the 
CAR1 (cyclic AMP receptor) gene |0|]; the expression 
of this gene leads to a membrane-localized fusion pro- 
tein which causes the cell to be fluorescently outlined, 
as shown. This new protocol for Dictyostelium develop- 
ment allows one to track cell motions in much greater 
detail than has been possible to date. 

It has been previously noted |l4) that coherent ro- 
tational motion can often be seen in three dimensional 
Dictyostelium mounds, albeit as a short-lived transient 
prior to cell-type sorting and tip formation at the mound 
apex. This motion has been attributed fill to cells mov- 
ing chemotactically to rotating waves of cyclic AMP. To 
test this hypothesis, we have repeated our experiments 
with a non-signaling strain of Dictyostelium |l6| . Aside 
from the need for a higher initial density to overcome the 
inability of the system to support long-range aggregation, 
the system behaves in a similar manner and produces 
rotating vortex structures. Thus, guidance via cAMP 
waves is not a necessary ingredient for organized rota- 
tional motion. Instead, we suspect a self-organization 
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FIG. 1. Pictures of the vortex state in the experiments. 
In (a) the entire cell body is fluorescent while in (b) only the 
cell membrane is fluorescent. The bar represents 100 [im in 
(a) and 10 fim in (b). 

of the system similar to that seen in molecular dynamics 
simulations of a confined set of particles undergoing dissi- 
pative collisions [|l7],[l8| , albeit with deformable particles 
and without an artificially imposed box. 

To verify that such a self-organized state is indeed pos- 
sible, we turn to a new model of cell motion. Following 
Glazier and Graner jl9|,|o| , we assume that cells move via 
(roughly) volume-preserving fluctuations of their shape. 
These fluctuations are driven by an effective free en- 
ergy pjj which incorporates specific biophysical mecha- 
nisms appropriate to Dictyostelium cells. First, cells stick 
to each other and tend to move in a manner which max- 
imizes cell-cell boundaries over cell-medium ones ]22j . In 
addition, each cell contains an active cytoskeleton which 
can generate forces by cycles of front protrusions and 
back retractions p3[. This force appears in our model 
as a potential which preferentially accepts fluctuations 
which move the cell in the direction of this propulsive 
force. To qualitatively verify this methodology, we have 
simulated the movement of a single self-propelled parti- 



FIG. 2. Movement of a single self-propelled cell (black) 
surrounded by non-propelled cells. Snapshots are taken every 
800 MCS. 

cle, as shown in Fig. ^. The sequence clearly illustrates 
the deformability of the cell and reproduces qualitatively 
the observed movement of a Dictyostelium cell. 

In the aforementioned previous work on flocking and 
related systems, a key insight is that local interactions 
can cause local velocity correlations which then causes a 
macroscopically ordered state to emerge. In our model, 
this suggests that one requires that a cell adjust its own 
propulsive force based on its interactions with its neigh- 
bors. We have investigated several mechanisms whereby 
this could be incorporated in our model. The most bio- 
logically plausible assumption is that the direction of the 
force is updated so as to better match the forces exerted 
by neighboring cells - essentially a minimal frustration 
hypothesis, which is in accord with what one observes 
directly in the video-microscopy. The results that fol- 
low were derived for this specific assumption, but depend 
very weakly on any details beyond the basic correlation 
idea. 

To investigate the possibility of coherent vortex mo- 
tion we started simulations with 100 square cells stacked 
in a square surrounded by medium. The force direction 
of each cell was chosen at random. A snapshot of a typi- 
cal final state is shown in Fig [| which shows the bound- 
aries of the cells and the force direction as a line which 
starts at the center-of-mass (CM) of the cell and which 
points in the direction of the force. The cells form a 
roughly circular patch and are rotating around the cen- 
ter of the patch. The final state is typically reached after 
a transient of 100-1000 Monte Carlo steps (MCS) @. 
Depending on the initial conditions, the cells will rotate 
either clockwise or anti-clockwise. Other than the sense 
of rotation and the duration of the transient phase, noth- 
ing depends on the starting state. Similar results were 
obtained for different numbers of cells. The existence of 
a localized rotating state is a robust consequence of the 
model, present as long as the cells are sufficiently cohe- 



FIG. 3. Snapshot of a typical final state of the model. 
The solid lines within each cell start at the CM of the cell 
and point in the direction of the force. The parameter values 
for this and all other figures in this Letter, are: N = 100, 
A target = 100, A = 10, Jcc = 5, J cm = 15, C = 1.0 and 
effective "temperature" T = 5. The lattice contains 200x200 
sites and the force direction is updated every 2.5 MCS. 

sive (Jem > Jcc)- Hence, we have proven that simple 
models based directly on the observed cell motions can 
indeed account for the self-localized vortex. 

To further test our model, we consider the angular ve- 
locity of cells as a function of radius. In the experiments, 
the angular velocity was obtained by direct tracking of 
cells. This tracking is greatly facilitated by the use of the 
cell-membrane outlining approach as shown in Fig. lb. 
In detail, we have taken six separate sequences of 15 min. 
each and measured the angular velocity every 8 s. During 
each sequence the radius of the individual cells changed 
little. Next, we grouped the data in radius intervals of 
4 fim and calculated the average velocity and the stan- 
dard deviation for each interval. The data is shown in 
Fig. ^ where the vertical bars represent one standard de- 
viation. The overall time scale corresponding to a MCS 
was adjusted to provide the best fit of the model (solid 
line) to the data; this yields 1 MCS as 0.006mm. As a 
consistency check, we note that this gives an isolated cell 
velocity of 8 ^m/min, which is very close to the exper- 
imentally reported value of 10 /im/min p4| . Although 
there is some way to go for a fully quantitative theory, 
our simple model does surprisingly well in capturing both 
the tendency of the cell-cell interaction to speed up the 
motion and the tendency of cells to slip as they move 
past each other, thereby limiting the angular velocity at 
large radii. 

What additional predictions does our model make? As 
the model proposes that cell-cell adhesion is the cause of 
the localized coherent state, mutant strains with reduced 
adhesion should not be able to form this structure. Also, 
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FIG. 4. The angular velocity as a function of the radius 
measured in the experiments (solid circles) and calculated us- 
ing the model (solid line). The overall time scale in the sim- 
ulations was adjusted to provide the best fit of the model to 
the data. 

mutants with weakened cytoskeletons would not be able 
to organize their motion. As chemical signaling is not 
necessary, disrupting the external cAMP concentration 
should have minimal effects. Finally, our model suggests 
that the time scale for organization should roughly be 
tens of minutes for aggregates with hundreds of cells; this 
is in qualitative agreement with our observations (data 
not shown). 

Vortex structures have been seen in other microor- 
ganism systems, namely the nutrient-limited spreading 
of the newly named bacterium Paenibacillus dendri- 
tiformis [ p5[ , and in Bacillus circulans p6] |. The mo- 
tion of bacteria occurs through flagella and is thus fun- 
damentally different from the amoeboid motion of Dic- 
tyostelium. However, it is tempting to speculate that 
in these cases as well velocity correlations induced by 
cell-cell interactions as well as cohesive forces may be 
enough to account for these structures [^7) . The fact that 
these disparate systems exhibit such strikingly similar 
non-equilibrium structures offers comfort to the physicist 
proposing a simplified model for an inscrutably complex 
biological process. 

In summary, we have documented the existence of a lo- 
calized, coherent vortex state in Dictyostelium. Further- 
more, we have argued both experimentally and via con- 
struction of a new model that the coherent motion is self- 
organized and not the result of a rotating chemical-wave 
guiding the motion. The advantage of this deformable- 
cell model is that it allows for the direct comparison of 
simulation with observation. In the future, we plan to ex- 
tend our calculations to the later stage processes of cell 
sorting and slug formation. 
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